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We develop a microscopic theory describing a quantum impurity whose rotational degree of free¬ 
dom is coupled to a many-particle bath. We approach the problem by introducing the concept of an 
“angulon” - a quantum rotor dressed by a quantum field - and reveal its quasiparticle properties us¬ 
ing a combination of variational and diagrammatic techniques. Our theory predicts renormalisation 
of the impurity rotational structure, such as observed in experiments with molecules in superfluid 
helium droplets, in terms of a rotational Lamb shift induced by the many-particle environment. 
Furthermore, we discover a rich many-body-induced fine structure, emerging in rotational spectra 
due to a redistribution of angular momentum within the quantum many-body system. 


The concepts of rotation and angular momentum are 
ubiquitous across quantum physics, whether one deals 
with the lifetimes of unstable nuclei [1], accuracy of 
atomic clocks [2], or electronic structure of defect cen¬ 
ters in solids [3]. Pioneered by the seminal works of 
Wigner [4] and Racah [5], the quantum theory of angular 
momentum evolved into a powerful machinery, commonly 
used to classify the states of isolated quantum systems 
and perturbations to their structure due to electromag¬ 
netic or crystalline fields [6, 7]. In “realistic” experiments, 
however, quantum systems are almost inevitably coupled 
to a many-particle environment and associated fields of 
excitations, which is capable of profoundly altering the 
physics of the system [8-10]. 

Although studying the effects of a fluctuating bath on 
the dynamics of quantum impurities represents a vast re¬ 
search field of its own [9, 10], an understanding of the 
many-body effects in the context of rotational degrees 
of freedom is still at its embryonic stage. As opposed 
to translational motion, whose coupling to a quantum 
field constitutes the well-studied “polaron problem” [11— 
14], quantum rotations in three-dimensional space are 
described by a non-Abelian SO( 3) algebra and possess 
a discrete spectrum of eigenenergies. This results in a 
tremendous complexity of the angular momentum prop¬ 
erties even for a few interacting particles [6]. On the 
other hand, unlike spin degrees of freedom [15], rota¬ 
tion is explicitly associated with an intrinsic motion of 
the system, whose coupling to a bath needs to be prop¬ 
erly accounted for by a microscopic theory. As a con¬ 
sequence, the description of rotating particles immersed 
into a many-body system cannot be reduced to any of 
the previously known impurity problems of condensed 
matter physics. 

Here we uncover rich physics associated with quan¬ 
tum rotation coupled to a many-particle environment. 
Using a combination of variational and diagrammatic 
techniques, we demonstrate that the problem can be de¬ 
scribed within the quasiparticle picture of an “angulon” - 
a quantum rotor dressed by a quantum many-body field. 
The angulon is a collective object, characterized by the 


total angular momentum of the system, of which it is an 
eigenstate. At the same time, due to the impurity-bath 
interactions, the angular momentum is shared within the 
many-particle system, which results in a peculiar behav¬ 
ior of the angular momentum eigenstates. In particular, 
we discover the rotational Lamb shift as well as a rich 
many-body-induced fine structure emerging due to the 
transfer of angular momentum between the impurity and 
the many-particle bath. The latter effect has no direct 
analogue in isolated atoms and molecules. We identify 
the regimes where the predicted effects are accessible by 
experiments on cold molecules in Bose-Einstein Conden¬ 
sates (BEC’s) and inside superfluid helium nanodroplets. 

As a first step, we derive the effective Hamiltonian 
of the angulon problem, starting from a microscopic de¬ 
scription. For the sake of clarity, we exemplify the deriva¬ 
tion by considering a linear-rotor molecule immersed into 
a weakly-interacting BEC. The resulting Hamiltonian, 
Eq. (2) is, however, applicable to a broad range of sys¬ 
tems consisting of a quantum rotor coupled to a many- 
body bath of harmonic oscillators, such as phonons in a 
liquid [16] or solid [17], as well as microwave photons or 
electromagnetic field noise [18]. 

The Hamiltonian for a molecular impurity in a homo¬ 
geneous BEC is given by H = Hbos T -Hi m p T H irn p-bos? 

where H bos = Ek e k«iUk + 9bb Ek,k', q 4'- q 4+q®k'ak 
describes a gas of bosons with a dispersion relation ek, 
whose contact interactions are given by the strength 
9bb > 0 [19]. In units where h= 1, the kinetic energy of 
a linear rotor impurity is given by = B J 2 , where B 
is the rotational constant and J is the angular momen¬ 
tum operator. Accordingly, the non-interacting impurity 
eigenstates, |j, m), are characterized by the angular mo¬ 
mentum, j , and its projection, m, onto the laboratory- 
frame z-axis, and form (2j + l)-fold degenerate multiplets 
with energies Ej = Bj(j + 1) [20, 21]. We consider a 
molecule whose linear motion is frozen, which is a good 
approximation for the experimentally relevant cases [16]. 

Anisotropic molecular geometry gives rise to 



anisotropic impurity-boson interactions, 

-^imp-bos = b^mp-bos(Q j 0)^k—q^k? (1) 

k,q 

which are explicitly dependent on the molecular orien¬ 
tation in the laboratory frame, as defined by the an¬ 
gles (0,0), see Fig. 1(a). Here the two-body inter¬ 
action, Ump-bos(q, 0 , 0 ) — ^F\R{0 0 ) Vimp-bos (**)] 5 is ob¬ 
tained by applying the rotation operator, R(6 ,0), to the 
molecular-frame potential, Ui mp _bos ( r/ 0 and T denotes 
the subsequent Fourier transform. For a linear rotor, 
the anisotropic potential can be expanded in spherical 
harmonics, bmp-bos(r') = Ea u \h( r ') Y \o(Q ',$')> where 
(r',©',^') are the boson coordinates in the molecular 
frame and u\ and f\(r') represent the strength and shape 
of the potential in the respective angular momentum 
channel A, cf. Fig. 1(b) [22]. 

After expanding the Hamiltonian in fluctuations 
around a homogeneous BEC of condensate density 
n [23-25], its bosonic part can be diagonalised us¬ 
ing the Bogoliubov transformation [19]. The result¬ 
ing theory is expressed in terms of bosonic quasipar¬ 
ticles (‘phonons’) with dispersion relation Wk- The 
corresponding creation and annihilation operators are 
conveniently expressed in the spherical basis, b\ Xfl = 

k(27r)~ 3 / 2 f dQk ^*^(0/0. Here k = |k|, and A 

and fi define, respectively, the phonon angular momen¬ 
tum and its projection onto the 2 :-axis. Up to a constant 
mean-field shift, the resulting effective Hamiltonian reads 


H — BJ + '5^Wkfys\ u pk\iJL 

k\(i 

+ Y, Ux(k) [y^( 9, 4>)b{ Xll + Y Xfl 0,4>)b kXfi 

kX/i 


( 2 ) 



FIG. 1. (a) The interaction of a quantum rotor with a 

quantum many-body system explicitly depends on the rotor 
angular coordinates, (0,0), in the laboratory frame, (b) The 
anisotropic rotor-boson interaction is defined in the rotor co¬ 
ordinate frame, r' = (r x , ©', $'). (c) Representation of Eq. (4) 
in terms of Feynman diagrams. 


On the other hand, the spherical harmonics appearing in 
Eq. (2) incorporate the non-Abelian algebra of angular 
momentum, which is not present in the Frohlich Hamil¬ 
tonian for a translationally moving impurity [14]. 

While exact solutions for Eq. (2) exist in the limit of 
a non-rotating molecule, this is no longer the case for fi¬ 
nite H, where the rotor becomes a dynamical degree of 
freedom. In order to get insight into the angulon proper¬ 
ties, including the non-perturbative regime, we introduce 
a variational ansatz for the many-body quantum state 
which is based on an expansion in bath excitations, 

IV>) = Z |0) I LM) + 1°) tf™>> (3) 

k\[i 

jm 


with = f dk and 

turn dependent interaction 

8 nk 2 tk 


U\ 


1 1/2 


the angular momen- 
given by U\(k) = 


f drr 2 f\(r)j\(kr) with j\(kr) the 


Wk (2A+1) 

spherical Bessel function [22]. 

It is important to note that, although we derived 
Eq. (2) in the context of a particular physical realisa¬ 
tion, an effective Hamiltonian of the same structure can 
be obtained for any quantum system whose angular mo¬ 
mentum is coupled to a bath of quantum oscillators. Hav¬ 
ing this in mind, we will hereafter regard Eq. (2) from a 
general perspective and study its generic properties both 
within and beyond the cold-atom regime. While the first 
two terms of Eq. (2) correspond to the bare kinematics 
of the impurity and the bath, the last term accounts for 
the absorption and emission of field quanta by the quan¬ 
tum rotor. As opposed to the spin-boson model [15], this 
term explicitly depends on the molecular angle operators, 
( 0 , 0 ), which is essential for the microscopic description of 
physical rotation and the emergence of angulon physics. 


where | 0 ) represents the vacuum of bath excitations and 
the angulon quasiparticle weight is given by the normal¬ 
isation condition, Z LM = 1 — J2k\j \Pk\j\ 2 . The many- 
body state (3) is an eigenstate of the total angular mo¬ 
mentum, L 2 |0) = L(L + 1)1-0), and its projection on the 
laboratory z-axis, L z |0) = M|0), which is incorporated 
by the Clebsch-Gordan coefficient C^f Xfl [ 6 , 7]. In the 
absence of external fields, the quantum number M is ir¬ 
relevant and will be omitted hereafter. 

Minimization of the energy, E = ( / ip\H |0)/(0|0), sub¬ 
ject to the constraint ( 010 ) = 1 , yields the variational 
energies, El, as solutions of the self-consistent equation, 
= 0 , where 

[Grm- 1 = [Gi(E )]- 1 - S l (E), (4) 

with [G^(E )]- 1 = —E + BL(L + 1) and 


£l(E) = E 

k\j 


2A + 1 U x(k) 2 

°L0,A0 

2 

47r Bj(j + 1] 

) — E + (dk 


(5) 
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FIG. 2. (a) The angulon spectral function, Al(E), as a function of the dimensionless density, n = n(mB )~ 3 ^ 2 , and energy 

E — E/B. The low- and high-density regimes can be realized with cold molecules trapped inside a BEC and superfluid 
helium droplets, respectively. The ground state, 0 q O , is stable; so is the state 1^ after it crosses the phonon threshold at zero 
energy, marked by the horizontal dashed line. Within the plotted range of n, the rest of the excited states have finite lifetimes, 
(b) Differential rotational Lamb shift for the lowest nonzero-L states, (c) Zoom-in illustrating the Many-Body-Induced Fine 
Structure (MBIFS) of the first kind, Ll,o —>• { L 7 0 , Lj 0 }, and of the second kind, L 7 0 —>> L^_ 1 1 . (d) Spectroscopic signatures 
of the MBIFS for the L — 1 state. The numbers indicate the corresponding values of ln[n]. Sharp features in panels (a), (c), 
and (d) were artificially broadened for better visibility. 


Eq. (4) is equivalent to the Dyson equation G ang = G° + 
G°EG ang for the interacting Green’s function, G^ ng (£ 1 ), 
with the self-energy, ^^(E), calculated from the dia¬ 
grammatic expansion shown in Fig. 1(c). As a conse¬ 
quence, our variational technique is equivalent to a dia¬ 
grammatic approach to the problem. As such, it allows 
to access not only the ground-state properties but also 
the entire excitation spectrum of the system by virtue of 
the retarded Green’s function, G r l l (E) = G^ g (E + i0 + ). 
The corresponding self-energy is evaluated in closed form: 

1 2 


Im£^ et 


(E) = y j e{E- Bj(j 


Aj/co 

2A 

x — 


1 )) 


rij 0 
°L0,A0 


-u x {k 0 y\{du k /dk) k = ko \~^ (6) 


with ReE^ t (£’) given by the Kramers-Kronig rela¬ 
tions [17]. Here fco gives the roots of E — + Bj(j + 1) = 

0. The theta-function represents the onset of phonon 
bands, shown in Fig. 2(d) (see below), which have been 
observed in experiments with helium nanodroplets [26]. 
The analytic results allow for the efficient calculation of 
the spectral function, Al(E) = Im[G^E 1 )], which en¬ 
compasses both the excitation spectrum as well as the 
quasiparticle properties of the angulon, and therefore is 
the central object to be studied below. 


Let us analyze the generic properties of Eq. (2), by 
considering a molecule immersed in a superfluid with 
dispersion relation c= ^/ck( e k + ^9bb n )^ where $bb = 
47rabb/^ [19]. Real atom-molecule potentials depend on 
the particular system and comprise high-order terms in 
their spherical harmonics expansion [27]. In Eq. (2), 
each term of the expansion leads to an effective phonon- 
mediated coupling between the bare rotational states. 
In order to keep the analysis of the Hamiltonian most 
transparent, we focus on the leading terms, A = 0,1, 
which suffices to reveal the intricate interplay between 
the anisotropic and isotropic interaction channels respon¬ 
sible for the emergent properties of the angulon. 

In Fig. 2(a) we show the resulting angulon spectral 
function, Al(E) = Al(E)B, as a function of the di¬ 
mensionless superfluid density, n = n{mB )~ 3 / 2 and en¬ 
ergy E = E/B. We choose Gaussian form factors, 
fx(r) = (27r) -3 / 2 e -r A 2r x ) ? an d interaction parameters 
uo = 1.75^i = 218H, ro = r\ = 1.5(mB) -1 / 2 in 
order to represent the strength and range of a typi¬ 
cal atom-molecule potential [27]. Furthermore we set 
<2bb = 3.3(mfi) -1 / 2 , which approximately reproduces 
the speed of sound in superfluid 4 He for a molecule with 
B = 2tt x 1 GHz [28]. 
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The long-lived angulon states correspond to the sharp 
features (dark shade) in Fig. 2(a), and we find that 
their energy decreases with growing density. How¬ 
ever, in addition to a uniform shift, also present in 
polarons [11-14, 29-32], in the angulon problem the 
anisotropic molecule-boson interaction leads to an addi¬ 
tional shift, whose magnitude depends on L. Such a shift 
has been observed in spectra of molecules embedded in 
superfluid helium nanodroplets [16, 33-35], however it 
has never been derived from a general microscopic the¬ 
ory. By analogy with quantum electrodynamics (QED) 
[36], this effect can be understood as “rotational Lamb 
shift” (RLS). Fig. 2(b) illustrates the RLS, defined as 
A^ LS = (El — Eo)/B — L(L + 1), for a few lowest angu¬ 
lar momentum states. The main feature is the change of 
sign of the RLS at the transition from the perturbative 
to the non-perturbative regime. As detailed below, this 
effect is closely related to the development of the many- 
body-induced fine structure of the second kind. We find 
that the RLS is quite sensitive to the anisotropy of the 
two-body potentials and therefore can be used as an ac¬ 
curate spectroscopic tool to measure the microscopic pa¬ 
rameters of the system. 

Even at low densities the magnitude of the RLS is suf¬ 
ficiently large to be measured in modern experiments on 
ultracold quantum gases. For example, for a ground- 
state molecule with B = 2tt x 1 GHz, trapped in a 23 Na 
BEC of density 10 15 cm -3 with abb = 60 a.u. [37] we ob¬ 
tain a Lamb shift on the order of 10 Hz. This prediction 
assumes molecule-atom interactions characterized by an 
s -wave scattering length a = —100 a.u., as well as an 
anisotropic A = 1 potential of comparable strength and 
finite range rq = 100 a.u. For ultracold molecules in 
highly-excited vibrational states with B = 2tt x 1 MHz 
we find an increase of the Lamb shift by an order of mag¬ 
nitude. 

In addition to the RLS, a rich fine structure of 
the angulon excitation spectrum unfolds with increas¬ 
ing interaction strength n. In spirit of molecular spec¬ 
troscopy [20], we classify the angulon states according to 
Z/yA, where in addition to the good quantum number, 
L, we introduced the “approximate quantum numbers”, j 
and A, indicating the dominant angular momentum con¬ 
tributions from the molecule and bosons, respectively. 
Starting from weak interactions, the quasi-free molecu¬ 
lar states, Ll, o, become dressed by bosonic excitations, 
which leads to a negative energy shift. While the ground 
angulon state, Oqo, is infinitely long lived, any of the 
L > 0 states has a finite lifetime due to the weak cou¬ 
pling to phonons with A > 0. For states with Zl ~ 1 the 
lifetime r is given by Eq. (6) as r _1 = ImE^ (El)- The 
^-function of Eq. (6) determines which decay channels 
are open at a given energy, while the Clebsch-Gordan 
coefficient gives the branching ratios between them. 

As interactions increase with h, two distinct features 
appear in the angulon spectrum which we refer to as 


Many-Body-Induced Fine Structure (MBIFS) of the first 
and second kind. The MBIFS of the first kind is related 
to a new angulon state emerging on top of the phonon 
continuum at positive energies when the interaction en¬ 
ergy suffices to probe the finite range of the isotropic part 
of the interaction, E ~ l/(2raro). As a result, the state 
splits into two angulon branches, Ll ,o — 
as shown in Fig. 2(c,d). An analysis of the self-energy in 
Eq. (5) shows that this effect is solely due to the isotropic 
part of the finite-range potential. The appearance of 
the state can be understood as a resonance in the 
many-body spectrum emerging due to coupling between 
the molecule and phonon states outside of the scatter¬ 
ing continuum. Within our model, the resulting excited 
state, L\ 0 , is stabilized since the scattering continuum is 
exponentially suppressed by the finite range interaction, 
ro (imposing a cutoff onto the bath spectral function, 
cf. the spin-boson model [15]). A similar stabilization 
mechanism appears e.g. for the repulsively bound atom 
pairs in optical lattices [38]. Furthermore, the state 
is reminiscent of the excited bound phonon-polaron state 
in the Holstein model [39, 40]. However, the formation of 
the upper-branch angulon L\ 0 competes with the cou¬ 
pling of the molecule to bosons with non-zero angular 
momentum. ‘On mass shell’ emission of such bosons out 
of the phonon vacuum, a process analogous to sponta¬ 
neous emission in QED, renders the L + state metastable 
and ultimately leads to its complete decoherence at large 
values of n. 

The phonon-mediated interactions between different 
rotor states, j, originate from anisotropies of the two- 
body potentials. Once these couplings reach a critical 
value, another peculiar effect occurs, which we refer to 
as MBIFS of the second kind. Here, strong entangle¬ 
ment with the bath leads to the disintegration of the 
lower-branch angulon, L^ 0 , into the phonon continuum 
belonging to the rotor state with j = L — 1. The L~ state 
then reappears at lower energy separated by a gap which 
scales with the anisotropic part of the potential, u\. An 
analysis of the occupation numbers f3^ X j of Eq. (3) re¬ 
veals that at the transition, approximately one quantum 
of angular momentum is transferred from the rotor to 
the many-body bath, L^_ 1 r Within the single¬ 

phonon excitation approach of Eq. (3), the reemerging 
angulon has an infinite lifetime as soon as it lies energet¬ 
ically below the scattering continuum of the L — 1 state 
as can be seen analytically from Eq. (6). The stability 
threshold for 1^ state is shown in Fig. 2(a) by the hor¬ 
izontal dashed line. Both the MBIFS of the first and 
second kind can be probed spectroscopically at the typi¬ 
cal parameters of molecules in helium nanodroplets [16]; 
their spectroscopic signatures are shown in Fig. 2(d). 

In summary, we developed a quasiparticle-based ap¬ 
proach to the redistribution of angular momentum be¬ 
tween a rotating impurity and a bosonic bath. Our the¬ 
ory sheds light on the properties of molecules trapped 
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in superfluid helium droplets from a general many-body 
perspective, and predicts a rich rotational fine-structure 
arising due to the interaction with a many-particle en¬ 
vironment. Accounting for the predicted effects might 
be key to the understanding of recent experimental data 
on non-adiabatic molecular excitations [41], which so far 
lacks even a qualitative explanation. In a broader per¬ 
spective, the angulon can serve as a building block for 
the description of few-, and many-molecule processes in 
the presence of an environment, such as reactivity and 
molecular collisions, and provide an efficient description 
of quantum rotor models [42] coupled to an external bath. 
The presented theory can be extended to account for 
complex molecular potentials and vibrations [27], exter¬ 
nal fields [21], the physics of rotons and maxons [28], 
multi-phonon excitations, as well as phonon interactions 
and decay [43, 44] which is expected to further enrich the 
observed phenomena. 
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Here we provide some technical details on the derivation of the angulon Hamiltonian, Eq. (2) of the main text. As 
described there, we start with the microscopic Hamiltonian describing the bosons, the impurity, and the impurity- 
boson interactions: H = Hbos + iAmp + ^imp-bos- The last term, given by Eq. (1) of the main text, explicitly 
depends on the orientation of the impurity in the laboratory frame, (0, <f >). As a first step, we expand the Hamiltonian 
in fluctuations around a homogeneous BEC of condensate density n and perform the Bogoliubov transformation, 
neglecting the term ~ b^b k, which is suppressed by 1 /y/n. In Cartesian coordinates, up to a constant mean-field shift 
~ V (0)n, the resulting Hamiltonian reads: 


H = BJ 2 + ^2 wfc&k&k + Timp- bos (k) @ •> 0) 

k k 



(1) 


where = f d 3 k/{2 tt) 3 . 


Impurity-boson interaction 

We now proceed with deriving the expression for Hi m p-bos(k, 0, </>)• As shown in Fig. 1(a) of the main text, in 
addition to the laboratory coordinate frame, (x,y,z), we define the molecular coordinate frame, (x',y',z'), needed 
to describe the microscopic molecule-boson interactions. In the molecular frame, an axially-symmetric interaction 
potential can be expanded over spherical harmonics as: 

bmp-bos(r') = X «a/a (r')Yx 0 (&, $')> (2) 

A 

where u\ and f\(r f ) represent, respectively, the strength and shape of the potential in the corresponding angular 
momentum channel A, cf. Fig. 1(b) in the main text, and (r', ©', 4> / ) represent the boson coordinates in the molecular 
frame. Next we transform Eq. (2) to the laboratory frame, (r, ©, <h), using Wigner rotation matrices [1, 2]: 

^0(0',$') = *) (3) 

u 

Here the Euler angles (0, 0 , 7 ) define the relative orientation of the two coordinate systems. For a linear molecule, the 
third angle, 7 , can be dropped. Then, using that 7>^j(</>, 0, 0) = yj 27+1 Ta m (^, (jf) [1, 2], we obtain: 

bmp-bo S (r, e,fr = j2 mxM h ( 4 ) 

X/i 

Here r = (r, 0, <h) gives the bosons’ positions in the laboratory frame. 

The interaction in the momentum space which appears in Eq. (1) is given by the Fourier transform of Eq. (4), 
evaluated using the expansion of e _2kr in terms of spherical harmonics: 

bmp-bos (k, f d 3 r bmp-bos(r, 0, 4>) e" ikr = ]T(2 7 T) 3/2 i~ x Vx(k)Y Xfl (e k , $ k )Y Xtl (§, <7 (5) 

J Xu 

where k = (fc, 0^,$^) is the momentum vector in the laboratory frame. This yields the expression V\(k) = 
u\2?! 2 /\/ 2 A + 1 f^° dr r 2 f\(r)j\(kr), with j\(kr ) the spherical Bessel function. 



Hamiltonian in the angular momentum representation 


It is convenient to express the creation and annihilation operators of Eq. (1) in the spherical basis using the following 
relations: 


&£ = ME £ i{ r a iv(0 fc , $*), (6) 


&U = (2^37 2 / sin0 fe % (7) 

Here the quantum numbers A and fi define, respectively, the angular momentum of the bosonic excitation and its 
projection onto the laboratory-frame z-axis. The transformations defined by Eqs. (6) and (7) correspond to the 
following commutation relations: 


[S k ,6| c ,] = (27r) 3 ( 5( 3 )(k-k') 


(8) 


— S(k k 


(9) 


After substituting Eq. (5) and (6) into Eq. (1) and integrating over the angles, the effective Hamiltonian in the 
spherical basis takes the following form: 


H = BY 


k\/i 


fjpkX^ 


E u *w f y A>’ 0) W 


k\(i 


( 10 ) 


where Yu = f dk and U\(k) = V\{k)ky/ne^Juo^ = u\ 


8nk 2 ek 
uk (2A+1) 


1/2 


/ drr 2 f x {r)jx{kr). 
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